IFUM-956-FT 



The impact of PDF and ol s uncertainties on 
Higgs Production in gluon— fusion at hadron colliders 



Federico Demartin a , Stefano Forte"' 6 , Elisa Mariani a , 
Juan Rojo 5 and Alessandro Vicini a ' 5 

a Dipartimento di Fisica, Universita di Milano, and 
b INFN, Sezione di Milano, 
Via Celoria 16, 1-20133 Milano, Italy 



Abstract 

We present a systematic study of uncertainties due to parton distributions and the 
strong coupling on the gluon-fusion production cross section of the Standard Model 
Higgs at the Tevatron and LHC colliders. We compare procedures and results when 
three recent sets of PDFs are used, CTEQ6.6, MSTW08 and NNPDF1.2, and we discuss 
specifically the way PDF and strong coupling uncertainties are combined. We find that 
results obtained from different PDF sets are in reasonable agreement if a common value 
of the strong coupling is adopted. We show that the addition in quadrature of PDF 
and a s uncertainties provides an adequate approximation to the full result with exact 
error propagation. We discuss a simple recipe to determine a conservative PDF+a s 
uncertainty from available global parton sets, and we use it to estimate this uncertainty 
on the given process to be about 10% at the Tevatron and 5% at the LHC for a light 
Higgs. 
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1 Introduction 



At the dawn of experimentation of LHC it is important to assess carefully the expected 
accuracy of standard candle signal and background measurements. Standard model Higgs 
production is clearly one such process. The main production mechanism for a scalar Higgs 
boson at the LHC is the gluon-fusion process (pp ->■ H + X) [I]. This process is also an 
especially interesting test case to study QCD uncertainties: on the one hand, it starts at 
0(a 2 s ) and it undergoes large O(al) corrections which almost double the cross section (2j 
[3j. On the other hand, it is driven by the gluon distribution, which is only determined 
starting at 0(a s ) (unlike quark distributions which can be determined from parton-model 
processes). Therefore, this process is quite sensitive both to parton distributions (PDFs) and 
a s uncertainties, and also on their interplay. 

Our knowledge of PDFs (see Refs. (4j|5)) and of a s (see Refs. (6j[T] ) have considerably 
improved in the last several years. However, they remain the main source of phenomenological 
uncertainty related to the treatment of the strong interaction: they limit the accuracy in a 
way which cannot be improved upon by increasing the theoretical accuracy. It is the purpose 
of this work to explore this uncertainty using Higgs production in gluon-fusion as a test case. 
Our goal is fivefold: 

• We would like to compare the procedure recommended by various groups to combine 
PDF uncertainties and a s uncertainties (and specifically Hessian-based approaches with 
Monte Carlo approaches) both in terms of procedure and in terms of results. 

• We would like to assess the impact of the correlation between the value of a s and the 
PDFs both when determining central values and uncertainty bands, and specifically 
understand how much results change when this correlation is taken into account in 
comparison to the case in which a s and PDF variations are done independently with 
results added in quadrature. 

• We would like to assess how much of the difference in results found when using different 
PDF sets (both for central values and uncertainty bands) is due to the PDFs, and how 
much is due to the choice of value of a s . 

• We would like to compare how much of the total uncertainty is due to a s and how much 
is due to PDFs. 

• We would like to arrive at an assessment of the value and PDF+a s uncertainty for this 
cross section and more in general at a procedure to estimate them. 

For each of these issues Higgs production through gluon-fusion is an interesting test case in 
that it is likely to provide a worst-case scenario: differences between results obtained using 
different PDF sets or following different procedures for the combination of uncertainties are 
likely to be smaller for many other relevant processes. For instance, in processes involving 
quark PDFs the correlation between PDFs and the value of a s is likely to be weaker, and 
thus the results found adding uncertainties in quadrature are likely to differ less from those 
obtained when the correlation between PDFs and a s is fully accounted for. 

The studies performed here will be done using PDF sets from the CTEQ, MSTW and 



NNPDF Collaborations, specifically the PDF sets CTEQ6.6 8 , MSTW08 |9| and NNPDF1.2 10 



11 . In order to account for the a s dependence, we will use PDF sets with varying a s which 
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have been published by CTEQ [§] and MSTW fl2], as well as NNPDF1.2 sets with vary- 



ing a s 13 which will be presented here for the first time. Comparison between CTEQ and 



MSTW on one side and NNPDF on the other side will enable us to contrast results obtained 
in the Hessian approach of the former with those found in the Monte Carlo approach of the 
latter. Computations will be performed at next-to-leading order (NLO) in the strong cou- 
pling a s , because, even though NNLO results for the process we study are available [14] , 
global parton fits with a full treatment of both DIS and hadronic data only exist at NLO (for 
instance, the MSTW08 set |:9|| only treats DIS fully at NNLO, while Drell-Yan is described 
using if -factors, and jets using NLO theory). There are of course several other sources of 
uncertainty on standard candles at colliders, such as electroweak uncertainties and further 
QCD uncertainties unrelated to PDFs, but their study goes beyond the scope of this work: 
here we concentrate on PDF uncertainties, and on the a s uncertainty which is tangled with 
them. 

The outline of the paper is as follows: in Sect. [2] we summarize the computation of the 
Higgs production cross section and the choice of value of a s . In Sect. [3] we discuss and compare 
PDF sets with varying a s , and specifically present the NNPDF1.2 sets with varying a s , which 
allow for a direct computation of the correlation between a s and the gluon. We then turn to 
a comparison of predictions obtained using different PDF sets: first, in Sect. [4] we study PDF 
uncertainties and compare predictions for the cross section and the PDF uncertainty on it 
obtained using different sets; then in Sect.[5]we discuss a s uncertainties and their combination 
with PDF uncertainties; finally in Sect. [6] we compare final results and discuss a procedure to 
construct a combined prediction from the available sets. Conclusions are drawn in Sect. [7j 
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2 The Higgs boson production cross section 



2.1 The hadron— level cross section 



The hadronic total cross section for the production of a Standard Model Higgs of mass m# 
via gluon-fusion at center-of-mass energy s/s is 

a(hi + h 2 -> H + X) = dxxdxz f aM (xx, f bM (x 2 , ii 2 F ) 



x I dz 5 ( z 



x x x 2 



<7ab(Z 



(1) 



where r# = m 2 H /s, hf is the factorization scale, f a ,hi(x , Hp) are the PDFs for parton a, (a = 
g, g, q) of hadron /ij, and o a b is the cross section for the partonic subprocess ab — > H + X at 
the center-of-mass energy s = X\X 2 s = m 2 H /z. The latter can be written as 



where the Born cross section is 



<r (0) (/4) zG ab (z,/j 2 R ) 



(2) 



(7 



(0) 



G> 2 (/4) 
512^^ 



(10 

7 



(3) 



and the sum runs over all quark flavors that appear in the amplitude Q^ ll \ with 



0(10 



In Eq. (|4]) we have defined 



m 



[2 -(!-%,) #(0,0,^)] 



'1 - % 9 



i/(0,0^) = -log 2 (^) 



(4) 



(5) 



with the standard notation for Harmonic Poly logarithms (HPLs). 
Up to NLO, 



G ab (z) = G ( ° b \z) + 



G 



(1), 
ab I 



(6) 



where a, 6 stand for any allowed parton. Exact analytic results, with the full dependence on 
the masses of the quarks running in the loop, have been obtained for the NLO coefficient 
function G^ b in 3] and more recently in terms of HPLs in [15] . 



All numerical 



results presented in this paper have been obtained using a code based on 



the expressions of Ref. 15 . We will consider only the gluon-gluon channel, and evaluate the 



total cross section at NLO-QCD, with the running of the strong coupling constant a s (fi 



R) 

implemented as discussed in Sect. 2J2 below. The default choice for the renormalization scale 
is [Mr = ran- The cross sections have been computed including the contributions due to the 
top and the bottom quark running in the fermion loop, with masses m t = 172 GeV and 
m b = 4.6 GeV; the value of the Fermi constant is G^ = 1.16637 10 -5 GeV~ 2 . The top mass 
has been renormalized in the on-shell scheme (3,15 
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2.2 The strong coupling a. 



, its most accurate 



Even though the strong coupling a s can be determined by a parton fit 12 
determination is arrived at by combining results from many high-energy processes, most of 
which do not depend on PDFs at all. A recent combined determination [7] is 

a s (M z ) = 0.1184 ±0.0007, (7) 

while the currently publisher^] PDG average [6] has a rather more conservative assessment of 
the uncertainty: 

a s (M z ) = 0.1176 ±0.002. (8) 

Both these uncertainties should be understood as one-cx, i.e. 68% confidence levels. 

Because it starts at order a 2 and it undergoes sizable O(af) corrections, the Higgs cross 
section is very sensitive to the central value of the strong coupling and it is thus important 
to compare results obtained using the same value of a s . However, different values of a s are 
adopted by various parton fitting groups. Specifically, for the PDF sets we are interested in, 
the reference values are 

a s (M z ) = 0.118 for CTEQ6.6, 

a s (M z ) = 0.119 for NNPDF1.2, (9) 
a s (M z ) = 0.12018 for MSTW08. 

Therefore, in order to obtain a meaningful comparison, we must study the dependence of 
results obtained with different sets when the value of a s is varied about the central values of 
Eq. Q. 

For comparison of relative uncertainties, which are less sensitive to central value of a s , 
but of course very sensitive to the range in which a s is varied, we will assume that the one-cx 
and 90% confidence level variations of a s are respectively given by 

A (68) a s = 0.0012 = 0.002/c 90 (10) 
A (9 V = 0.002, (11) 

where ego = 1.64485... is the number of standard deviations for a gaussian distribution that 



correspond to a 90% C.L. interval. Our choice Eqs. (10 11) is thus intermediate between the 
choices in Eqs. ([8]-[7]). A reassessment of the value a s and and its uncertainty go beyond the 
scope of this work: these values are chosen as a reasonable reference, and ensure that our 
results will also be valid for other reasonable choices in the same ballpark. 

An important subtlety is the number of active flavours in the running of the strong cou- 
pling, as implemented by various PDF analyses. Indeed, QCD calculations are usually per- 



formed in a decoupling scheme 17 , in which heavy flavours decouple at scales much lower 
than their mass. When studying a process like Higgs production for a wide range of the Higgs 
mass one must then specify what to do above top threshold, and specifically fix the scale at 
which heavy quarks decouple, which amounts to a choice of renormalization scheme. This 
choice should be used consistently in the running of a s , the evolution equations for PDFs, and 
the computation of hard matrix elements. In Refs. [3 15 , a scheme in which the number of 
active flavours becomes rir = 6 at Q 2 = m 2 is adopted; this scheme is also used by NNPDF. 



1 The 2009 web PDG update 16 no longer provides a combined determination of a s , and refers to Ref. [7]. 
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Figure 1: Ratio of the strong coupling determined with a number of flavours that varies from 
rif = 5 to rif = 6 at Q > nit to that with rif = 5 at all scales. 



However, CTEQ and MSTW instead use a scheme in which rif = 5 even when Q 2 > mf both 
in the running of the strong coupling and in PDF evolution. 

The effect of this on the running of a s is small but non-negligible: in Fig. [T] we plot 
the ratio of the variable-flavor a s to the fixed-flavor one, and show that for Q = 2m t the 
discrepancy is already of order of 1%, and thus the effect on a quantity which depends on 
higher powers of a s accordingly larger. Of course, this scheme dependence cancels to a large 
extent once PDFs and hard cross sections are consistently combined, and a fully consistent 
comparison would require the use of the same scheme everywhere. This is difficult in practice 
because of the different choices adopted by NNPDF on the one hand, and CTEQ and MSTW 
on the other hand. For the sake of comparisons below, we will use the NLO Higgs cross 
section from Refs. [3jfl5] and, consistently a s which runs with rif that varies at Q = m t : this 
is then also consistent with the evolution equations used to construct the NNPDF set, but 
not with those used to construct the CTEQ and MSTW set. It should be born in mind that 
this incomplete cancellation of the scheme dependence above the top threshold may lead to 
a spurious difference between central values at the percent level between NNPDF and other 
groups. 
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Figure 2: The ratio of the central gluons obtained in NNPDF1.2 fits when a s is varied, divided by 
the reference NNPDF1.2 gluon at the initial evolution scale Qq = 2 GeV 2 . The comparison is shown 
both in a linear (left) and logarithmic (right) scales. One-u uncertainty bands are shown for the 
central, highest and lowest values of a s . 



3 Parton sets with variable strong coupling 
3.1 PDFs and a s in a Hessian approach 



PDF sets with varying a s have been presented by the CTEQ [8] and MSTW 12 collabo- 
rations. Specifically, CTEQ has released the CTEQ6.6alphas sets, which add to the central 
CTEQ6.6 fit H|, which has a s = 0.118, four more sets with a s = 0.116, 0.117, 0.119, 0.120. 
In these fits, a s is taken as an external parameter which is fixed each time to the given value 
along with all other physical parameters in the fit (such as, say, the fine-structure constant or 
the W mass). However, eigenvector PDF sets for the computation of PDF uncertainties are 
only provided for the central CTEQ6.6 set. Therefore, it is possible to study the correlation 
between the value of a s and PDFs, but not the correlation with their uncertainties. 

The MSTW collaboration instead has performed a simultaneous determination of PDFs 
and a s , which is thus not treated as a fixed external parameter, but rather as a fit parameter, 
which leads to the central value quoted in Eq. Furthermore, MSTW has also released sets 
of PDFs, analogous to the CTEQ sets discussed above, in which a s is taken as an external 
parameter, and varied in steps of 0.001 for 0.110 < a s < 0.130. The sets in which PDFs 
and a s are determined simultaneously may be used for a determination of the correlation 
between the value of a s and both PDF central values and PDF uncertainties, though with 
the limitation that the value and uncertainty on a s found in the fit must be used. 

3.2 PDFs and a s in a Monte Carlo approach 

In the NNPDF parton determination, a s is taken as a fixed parameter as in the CTEQ fits; 
the correlation between NNPDF1.2 PDFs and the value of a s has been discussed recently 



in Ref. 18 . For the present work, we have constructed a family of NNPDF1.2 PDF sets 
using different values of a s . The very recent NNPDF2.0 PDF set [19] also includes PDFs 
determined with different fixed valued of a s ^\ 

2 The NNPDF1.2 sets with variable a s are available upon request. The NNPDF2.0 sets with variable a s 
are available from the webpage of the NNPDF Collaboration, htpp://sophia. ecm.ub.es/nnpdf, and will 
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We have repeated the NNPDF1.2 PDF determination with a s varied in the range 0.113 < 
a s < 0.128 and all other aspects of the parton determination unchanged: for each value of a s 
we have produced a set of 100 PDF replicas. In Fig. [2] we show the ratios of the central gluons 
obtained in these fits compared to the reference NNPDF1.2 gluon with a s (Mg) = 0.119, 
together with the PDF uncertainty band which corresponds to the reference value. 

The qualitative behaviour of the gluon in Fig.|2]can be understood as follows. In NNPDF1.2, 
the gluon is determined by scaling violations of deep-inelastic structure functions, i.e. mostly 
from medium and small x HERA data, with the large x gluon constrained by the momentum 
sum rule. With a given amount of scale dependence seen in the data, smaller values of a s 
require a larger small x gluon, and thus because of the sum rule a smaller large x gluon. 
In global fits [8, 12, 19 the behaviour is essentially the same, up to the fact that some extra 
constraint on the large x gluon is provided by Tevatron jet data, as quantified in [19] . 

The size of this correlation of the gluon with the value of a s shown in Fig. [2] is clearly 



statistically significant; however, when a s is varied within its uncertainty range, Eq. (10), the 



change in gluon distribution is generally smaller than the uncertainty on the gluon itself. It 
is interesting to note that the size of the uncertainty for values of a s which are away from the 
best fit is often larger than the uncertainty when a s is at or close to its best fit value: this 
is to be contrasted to what happens in a Hessian approach, where linear error propagation 
inevitably implies that the PDF uncertainty shrinks when a s moves away from its best-fit 



value 12 



Given sets of replicas determined with different values of a s , it is possible to perform 
statistics in which a s is varied, by assuming a distribution of values for a s . For instance, the 
average over Monte Carlo replicas of a general quantity which depends on both a s and the 
PDFs, J 7 (PDF, a s ) can be computed as 



rep N 



i yet ^"jj 



J) fJj) 



(12) 



where PDF^^ stands for the kj- th replica of the PDF fit with a s = a$\ and the numbers 



of replicas for each value of a s in the total sample are determined by the probability 



distribution of values of a s , with the constraint 



Ma 



N, 



rep 



U) 

s 

rep ' 



(13) 



3=1 



Specifically, assuming that global fit values of a s (such as Eqs. (|7j(8])) are gaussianly dis- 
tributed, the number of replicas is 



Nrep OCexp 



\ 



(14) 



with the normalization condition Eq. (13) 



also be available through the LHAPDF interface. 
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Figure 3: Comparison of the gluon PDFs from MSTW08 (top), CTEQ6.6 (center) and NNPDF1.2 
(bottom) at the scales Q 2 = 4 GeV 2 and Q 2 = 10 4 GeV 2 as a s is varied, normalized to the 
corresponding central sets, determined with the value of a s listed in Eq. Q. The one-a uncertainty 
band for the central set is also shown in each case. 
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Figure 4: Comparison of the gluon PDFs from MSTW08, CTEQ6.6 and NNPDF1.2 at the scales 
Q 2 = 4 GeV 2 (upper plots) and Q 2 = 10 4 GeV 2 (lower plots), each determined with the value of a s 
listed in Eq. ([9j), all normalized to the central MSTW08 set. 



3.3 Comparison of global PDF sets with variable a s 

The dependence of the gluon distribution on the value of a s is summarized in Fig. [3j where 
results obtained using the CTEQ6.6, MSTW08 and NNPDF1.2 ot s series are compared, both 
at the scale Q 2 = 4 GeV 2 , close to the scale at which PDFs are parametrized, and at the high 
scale Q 2 = 10 4 GeV 2 , typical of electroweak final states. The three central sets of Fig. [3] are 
then compared directly in Fig. |4j The gluon luminosities computed from these three central 
sets are then finally compared in Fig. [5] at Tevatron and LHC energies. The luminosities 
are plotted as a function of mj for given energy using the leading-order kinematic relation 



x = m 2 H /s. It should be born in mind that, because of soft-gluon dominance 20 , the NLO 
contribution is strongly peaked at the endpoint, and thus probes mostly the luminosity at the 
same value of x as the LO. 

Notable features of these comparisons are the following: 



The same correlation pattern between the gluon and a s discussed for NNPDF in Sect. 
is also apparent for other sets. 

At Q 2 = 4 GeV 2 uncertainties at very small x for MSTW and NNPDF are large enough 
to swamp the a s dependence. This does not happen for CTEQ, likely due to the more 
restrictive gluon parametrization used. However, at Q 2 = 10 4 GeV 2 the region at which 
uncertainties blow up is pushed to much smaller values of x. 
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Figure 5: Comparison of the gluon luminosities from MSTW08, CTEQ6.6 and NNPDF1.2 at 
the Tevatron (upper left plot), LHC 7 TeV (upper right plot) and LHC 14 TeV (lower plot) and 
Q 2 = 10 4 GeV 2 , each determined with the value of a s listed in Eq. all normalized to the central 
MSTW08 set 



• Uncertainties at large x also swamp the a s dependence, due to the scarcity or (for 
NNPDF1.2) lack of data in this region. 

• Even in the medium x region where the gluon is best known PDF uncertainty bands are 
rather larger than the gluon variation due to the variation of a s within its uncertainty 
Eq. or even Eq. (§). 

• The three gluons in Fig. [4] overlap to one-a in most of the kinematic region. However, at 
low scale they disagree significantly at large x, and at high scale they disagree, though 
by less than about two a, both at very large and very small x. 

• Once gluons are convoluted into a parton luminosity, most of the disagreements seen in 
Fig. [4] are washed out: indeed, the parton luminosities Fig. [5] computed from CTEQ6.6, 
MSTW08 and NNPDF1.2 all agree within uncertainties, in the sense that their one-cr 
error bands always overlap (though sometimes just about). 
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NNPDF1 .2 - Correlation of g(x,Q 2 ) and a(M z 2 ) 



NNPDF1 .2 - Correlation of PDFs(x,Q 2 ) and a(M z 2 ) 
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Figure 6: The correlation Eq. (15) between PDFs and a s (Mz) as a function of x. Left: gluon PDF 



at Q 2 = 2 GeV 2 and Q 2 = 10 4 GeV 2 ; right: singlet, triplet and total valence PDFs for Q 2 = 2 GeV 2 . 



3.4 Correlation between PDFs and 



a. 



Correlations between different PDFs, or between PDFs and physical observables, have been 
computed by CTEQ using a Hessian approach [8~], and by NNPDF using a Monte Carlo 
approach 
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Within a Monte Carlo approach it is in fact easy to estimate the correlation 
between any pair of quantities by computing their covariance over the Monte Carlo sample. 
Statistics involving the value of a s can then be performed provided only replicas with different 



values of a s are available, as discussed in Sect. 3.2 above, Eq. (12). 



Indeed, in a Monte Carlo approach the correlation between the strong coupling and the 
gluon (or any other PDF) is given by 



p[a s (M 2 ),g(x,Q 2 )] 



a s (Ml) g (x, Q 2 )) rcp - (a. (Mf)> rcD (g (x, Q 2 )) 



rep 



rep 



a 



as (Ml) a 9{x,Q 2 ) 



(15) 



for each value of a s according to Eqs. (13 14) above. 



CO 



where the distribution of values of a s is automatically reproduced if one picks N™^ of replicas 



Our results for the correlation coefficient between the gluon and a s (Mz) as a function of 
x, computed using Eq. (|l5|), with the NNPDF1.2 PDFs of Sect. |3^2|and Eq. Q for a s , both 
at a low scale Q 2 = 2 GeV 2 and at a typical LHC scale Q 2 = 10 4 GeV 2 , are shown in Fig. [6j 
It is interesting to note how evolution decorrelates the gluon from the strong coupling. We 
also show in Fig. [6] the correlation coefficient for other PDFs: as expected for the triplet and 
valence PDFs it is essentially zero, that is, in NNPDF1.2 these PDFs show essentially no 
sensitivity to a s . The correlation coefficient Fig. [6] quantifies the qualitative observations of 
Fig.|2} 
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Figure 7: Cross sections for Higgs production from gluon-fusion the Tevatron (top), 
LHC 7 TeV (center) and LHC 14 TeV (bottom). All uncertainty bands are one-o" PDF 
uncertainties, with a s fixed according to Eq. Q. The left column shows absolute results, 
the central column results normalized to the MSTW08 result, and the right column results 
normalized to each group's central result. 



4 The cross section: PDF uncertainties 

We now turn to the computation of the Higgs production cross section and associate un- 
certainty band due to PDF variation, with a s kept fixed at each group's preferred value, 
as given by Eq. ([9]), and uncertainties consistently determined as one-cx intervals using each 
groups recommended method, hence in particular Hessian methods for CTEQ and MSTW 
and Monte Carlo methods for NNPDF. We will not address here the issue of comparison of 
the various methods, and we will simply take each group's results at face value, in particular 
by taking as a 68% confidence level the interval determined as such by each group. When 



comparing results it should be born in mind that, as discussed in Sect. 3.3 the cross sections 
probe the gluon luminosity at an essentially fixed value of x = m 2 H /s. 

4.1 Comparison of cross sections and uncertainties 



In Fig. [7] we compare the cross sections at Tevatron and LHC (7 TeV and 14 TeV) energies as 
a function of the Higgs mass; each cross section is calculated using the best PDF set of each 
group and the corresponding value of a s of Eq. ^ in the determination of the hard cross 
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section. Central values can differ by a sizable amount. Discrepancies are due to three distinct 
reasons: the fact that the hard cross sections (independent of the PDF) are different because 
of the different values of a s ; the fact that the PDFs are different because they depend on 
a s as shown in Figs. |3j lastly the fact that even when the same value of a s is adopted PDF 
determination from different groups do not coincide. 

As discussed in Sect. [3j the first two effects tend to compensate each other at small x 
because of the anticorrelation between the gluon and a s (see Fig. [6|, while at large x they go 



in the same direction. As we shall see in Sect. 5.2 the transition between anticorrelation to 



correlation happens for LHC 7 TeV for intermediate values of the Higgs mass. 

The relative impact the first effect (which affects the hard cross section) and of the second 
two combined (which affect the parton luminosity) can be assessed by comparing the cross 
sections of Fig. [7] with the luminosities of Fig. |5j for instance, for m# = 150 GeV at the 
LHC (7 TeV), the MSTW08 cross section is seen in Fig. [7] to be by about 7% higher than 
the CTEQ6.6. Of this, Fig. [5] shows that about 3% is due to the different parton luminosity, 
hence about 4% must be due to the choice of a, in the hard matrix element. In Sect. 5.2 we 



will determine this variation directly (Fig. 10) and see that this is indeed the case. 

Because the cross section starts at order or s (mH)i with a NLO iT-factor of order one, we 
expect a percentage change Aa s in a s to change the cross section by about 2.5Aa s , which 
indeed suggests a 4% change of the hard cross section when a s is changed from the M STW08 
to the CTEQ6.6 value. In fact, comparison of Fig. [7]to Fig. [5] shows that this simple estimate 
works generally quite well: the difference in hard matrix elements is 2.5Aa s so 4% when 
moving from the MSTW08 to the CTEQ6.6 value, with the rest of the differences seen in the 
cross sections in Fig. [7] due to the gluon luminosities displayed in Fig. [5j Because the latter 
are compatible to one-cr, this is already sufficient to show that nominal uncertainties on PDF 
sets are sufficient to accommodate the different central values of NNPDF1.2, CTEQ6.6 and 
MSTW08 once a common value of ot s is adopted. 

Uncertainties turn out to be quite similar for all groups and of order of a few percent, 
with uncertainties largest at the Tevatron for large Higgs mass. The growth of uncertainties 
as the Higgs mass is raised or the energy is lowered is due to the fact that larger x values are 
then probed, where knowledge of the gluon is less accurate, as shown in Fig. |4| 

4.2 The uncertainty of the PDF uncertainty bands 

In order to answer the question of the compatibility of different determinations of PDFs or of 
physical observables extracted from them it is important that the uncertainties are provided 
on the quantities which are being compared. Whereas this is standard for central values, it 
is less frequently done for uncertainties themselves. A systematic way of doing so in a Monte 
Carlo approach has been introduced in Refs. 10 19] (see in particular Appendix A of [19] ^ 



the difference between two determinations of a central value or an uncertainty is compared 
to the sum in quadrature of the uncertainty on each of the two quantities. Only when this 
ratio is significantly larger than one is the difference significant. 

It is important to observe that when addressing the compatibility of two determinations 
two inequivalent questions can be asked: whether the two determination come from sta- 
tistically indistinguishable underlying distributions, or whether they come from statistically 
distinguishable distributions, but are nevertheless compatible. To elucidate the difference, 
consider two different determinations of the same quantities based on two independent sets 
of data. If the data sets are compatible, the two determinations will be consistent within 
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Figure 8: The uncertainty on the Higgs cross section determined using NNPDF1.2 with different 
values of a s . The width of the band shown denotes the statistical uncertainty on the mean uncer- 
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tainty of the sample, c[c 2 ], Eq. (17). Note that the uncertainty on the uncertainty of the PDFs is 
Y^/Vrep = 10 times larger. 



uncertainties in the sense that the central values are compatible within uncertainties (i.e. the 
central values will be distributed compatibly to the given uncertainty if many determinations 
are compared). However, the underlying distribution of results will in general be different: for 
example, if one of the two datasets is more accurate than the other, the two distributions will 
certainly be statistically inequivalent (they will have different width), yet they may well be 
compatible. We expect complete statistical equivalence of two different PDF determinations 
only if they are based on the same data and methodology: for example, NNPDF verifies that 
its PDF determination is statistically independent of the architecture of the neural networks 
which are used in the analysis. However, when comparing determinations which use differ- 
ent datasets, PDF parametrization, minimization algorithm etc. we only expect them to be 
compatibile, but not statistically indistinguishable. 

The determination of the of the uncertainty on the uncertainty appears to be nontrivial 
in a Hessian approach, and it has not been addressed so far in the context of Hessian PDF 
determination to the best of our knowledge. We will thus only discuss it within a Monte 
Carlo framework. It seems plausible then to assume that uncertainties on uncertainties are 
of similar size in existing parton fits, so that a difference in uncertainties between two given 
fits is only significant if it is rather larger than about a/2<t[<7 2 ], where cr[<7 2 ] is the typical 
uncertainty which we will now determine in the Monte Carlo framework using NNPDF1.2 
PDFs. 

In a Monte Carlo approach, the uncertainty which should be used when assessing statisti- 
cal equivalence of two quantities determined from two sets of A^ rep replicas is the uncertainty of 
the mean of the replica sample, while the uncertainty to be used when assessing compatibility 
is the uncertainty of the sample itself, which is larger than the former by a factor iV r ep |10p9] . 
The former vanishes in the limit N Iep , while the latter does not: indeed, if the two replica sets 
come from the same underlying distributions all quantities computed from them should coin- 
cide in the infinite-sample size limit, while if they merely come from compatible distributions 
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even for very large sample quantities computed from them should remain different. 
Given a set of iV re p replicas, the variance of the variance can be determined as [6] 



(7 [a 



rep 



m A [q\ 



N rcp — 3 ^_2^2 



rep 



(16) 



where a 2 and WI4 are respectively the variance and fourth moment of the replica sample. 
Equivalently, we may determine c 2 [cr 2 ] by the jackknife method (see e.g. Ref. |21 



i.e. 



removing one of the replicas from the sample and determining the iV re p variances 



»?(*) = jf 



rep 



/ii(x)) 



2. 



I TV 

rep • 



(17) 



The variance of the variance is then given by 



iVr, 



E W-°t 



where 07 is given by Eq. (17) and 



(19) 



rep 



The quantity o\a 2 \ estimated using Eq. (16) or equivalently Eq. (18) provides the uncer- 



tainty of the sample, and indeed it vanishes in the limit iV r ep — > 00: it should thus be used 
to assess statistical equivalence. The quantity which should be used to assess consistency is 
JN rep a[a 2 



In Fig. |8|we plot cr[a 2 ] determined using Eq. (18) using the NNPDF1.2 sets with different 
values of a s . We can immediately conclude from this figure that the uncertainties on the 
cross sections does not change in a statistically significant way when a s is varied within its 
uncertainty Eq. ^ or even Eq. Q: a s has to be be varied by about three times the uncertainty 
Eq. (Jsj> for the uncertainty on the cross section to change in a statistically significant way. 
Incidentally, this plot also shows that uncertainties are not necessarily smaller (and in fact 



are mostly larger) when a s is varied away from its preferred value, as already seen in Sect. 3.2 
and Fig. |2j 



When assessing compatibilit 
certainty bands shown in Fig. 



(as opposed to statistical equivalence), the size of the un- 

This rescaled 



8 must be rescaled by a factor y/N^ 



rep 



10. 



uncertainties on the uncertainty are then more then half of the uncertainty itself. This means 
that all uncertainties shown in Fig. [7] for the three PDF fits under investigation are in fact 
compatible with each other at the one-a level, since they differ at most by a factor two 
(NNPDF vs. MSTW at LHC 14 TeV and the lowest values of m#), and in fact usually rather 
less than that. 
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5 The cross section: a s uncertainties 



We have seen in the previous Section (see Fig. [7]) that once differences in hard matrix elements 
due to the different choice of a s are accounted for, cross sections computed using different 
parton sets agree to one-a because parton luminosities do. However, parton luminosities 
compared in Fig. [5] were determined using the different values of a s Eq. ^ . 

For a fully consistent comparison, we must determine central parton luminosities (and 
thus cross sections) with a common value of a s , thereby isolating the differences which are 
genuinely due to PDFs. The uncertainty related to the choice of a s must be then determined 
by variation around the central value. This then raises the question of the correlation between 
this a s variation and the values of the PDFs and their uncertainties. We now address all these 
issues in turn. 



5.1 The cross section with a common a s value 



The cross sections obtained from the three PDF sets under investigation at the same value of 
a s are compared in Fig. [9j where the ratios of cross sections computed using in the numer- 
ator and denominator two different PDF sets (NNPDF1.2 vs. MSTW08 and CTEQ6.6 vs. 
MSTW08) are shown for different choices of a s . 

Comparison of these ratios to the PDF uncertainties shown in Fig. [7] show that they are 
typically of a similar size: namely, two-three percent at the LHC, while at the Tevatron 
they are about twice as large for light Higgs mass, rapidly growing more or less linearly, 
up to around 10% around the top threshold. The largest discrepancy in comparison to the 
PDF uncertainties is found at the lowest Higgs masses at the Tevatron between the CTEQ 
and MSTW sets, whose uncertainty bands barely overlap there. This can be traced to the 
behaviour of the gluon luminosities of Fig. [5] 

Hence, the spread of central values obtained using the PDF sets under investigation is 
consistent with their uncertainties, which are thus unlikely to be incorrectly estimated. Of 
course, it should be born in mind that these uncertainties are in turn estimated with a finite 



accuracy, unlikely to be much better than 50%, as seen in Sect. 4.2 



5.2 Uncertainty due to the choice of a s 

The simplest way to estimate the uncertainty due to a s is to take it as uncorrelated to the 
PDF uncertainty, and determine the variation of the cross section as a s is varied. This, in 
turn, can be done either by simply keeping the PDFs fixed, or else by also taking into account 
their correlation to the value of a s discussed in Sect. |3j i.e. by using for each value of a s the 
corresponding best-fit PDF set. In either case, the uncertainty on the cross section due to 
the variation of a s is 

{Aa)i s = a{a s ±Aa s )-a{a s ), (20) 

with, in the two cases, the cross section computed either with a fixed set of PDF, or with the 
sets of PDFs corresponding to the three values a s ± Aa s . 



We have determined (Acr)^ s Eq. (20) using the central value of a s of Eq. (10) and Aa s = 
0.002, which would correspond to a 90% C.L. variation of a s according to Eq. (11), and it 
is (almost exactly) equal to the difference between the preferred values of a s for CTEQ6.6 
and MSTW08, Eq. Q. For NNPDF1.2 a PDF set with a s = 0.117 is not available, hence 
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Figure 9: Ratio of the Higgs production cross sections determined using PDFs from different groups, 
but obtained with the same value of a s : CTEQ6.6/MSTW08 (left) and NNPDF1.2/MSTW08 
(right). Results are shown for different values of a s , and for the Tevatron (top), LHC 7 TeV 
(center) and LHC 14 TeV (bottom). 
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Figure 10: Relative variation Aa/a Eq. (20) of the cross section due a variation Aa s = 0.002 
about the preferred value Eq. (10) adopted by each PDF set. Results are shown both with 
PDFs kept fixed (dashed bands), or with the best-fit PDF set corresponding to each value of 
a s (solid bands). 



the lower cross section in the case in which the PDFs are varied has been determined by a 
suitable rescaling of that which corresponds to a s = 0.116. 



Results for LHC 7 TeV are shown in Fig. 10 upper and lower cross sections correspond 
to the upper and lower values of a s , but the variation turns out to be symmetric about the 
central value to good approximation. If the PDF is kept fixed, we find (Aa)^Jcr « 1.04 when 
Aa s = 0.002, i.e. a variation of about 4%, in excellent agreement with the simple estimate 



of Sect. 4A If the PDFs are also varied, the width of the uncertainty band on the cross 
section for light Higgs becomes smaller, because of the anticorrelation of the small x gluon to 
a s discussed in Sect. [3] (see Fig. pi), but it becomes wider as the the Higgs mass is raised and 
larger x values are probed, for which the correlation of gluon to a s is positive. 



5.3 Impact of a s on uncertainties 

So far we have considered the uncertainties due to PDFs and due to the value of a s separately. 
However, as seen in Sect. [3j when a s changes, not only the central values but also the uncer- 
tainties on PDF change, and this leads to a correlation between PDF and a s uncertainties. 
The effects of this correlation on the determination of combined PDF+a s uncertainties is 
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Figure 11: The combined PDF+a s relative one-cr uncertainty on the Higgs cross section with 
NNPDF1.2 (left) and MSTW08 (right). The three bands correspond to results obtained keeping into 
account the full correlation between a s and PDF uncertainties (exact), by adding in quadrature PDF 
uncertainties and a s uncertainties in turn determined by keeping into account the correlation between 
a s and PDF central values (quadrature), and finally by adding in quadrature PDF uncertainties and 
a s uncertainties determined with PDFs fixed at their central value (fixed PDF+Aa s ). The central 
values of a s are given by Eq. Q, and its one-cr range is in Eq. (10) in all cases except the MSTW08 
exact for which it is given in Eq. (22). 



likely to be moderate: it is a higher order effect, and the dependence of uncertainties on a s 
is weak, especially if compared to the their own uncertainty, see Fig. |8j 

We will now quantify this correlation by computing the total PDF+a s uncertainty when 
the correlation is kept into account, and comparing results to those obtained when PDF and 
a s uncertainties are added in quadrature. As discussed in Sect. |3j the correlated uncertainty 
can be determined both in the Hessian approach using MSTW08 and in a Monte Carlo 
approach using NNPDF. In the MSTW methodology this is done by relying on a simultaneous 
determination of PDFs and a s . As a consequence, the value and range of variation of a s must 
be those obtained in this determination. This will not hamper our analysis in that the 
MSTW08 value and range for a s of Ref. [12] are close to those under consideration. With 
the NNPDF methodology we are free to choose any value and range for a s , inasmuch as the 
corresponding Monte Carlo PDF replicas are available. 

We have thus computed joint PDFs+a s uncertainties on the Higgs cross section. For 



MSTW08 we have followed the procedure of Ref. 12 : the total upper and lower (generally 



asymmetric) uncertainties on an observable F are determined as 

(AF) PDF +Qs = m ^( {F «s {So) + {AF « h) + }) _ F a° {So) (21) 

(af) pdf +Qs = F *° {So) _ min ({F ^( 5o ) _ (AF p a £, F )_}, ) 

where F aa (So) is the observable computed using the central PDF set So an d the value a s 
of the strong coupling, (AFpf )F )± is the PDF uncertainty on the observable for given fixed 
value of a s , as determined from the Hessian PDF eigenvectors [9,12 , and the maximum and 



minimum are determined from a set of five results, each computed with one distinct value of 
a s (central,± half confidence level, ± confidence level). The value of a° s is as given in Eq. Q, 
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while its one-cx upper and lower variations are 



A(68) _+0.0012 
^ "a — -0.0015 



(22) 



For NNPDF, the uncertainty is simply given by the standard deviation of the joint distri- 
bution of PDF replicas and a s values 



AF PDF+o 



a F 



N a N IC S P 



rep j=l fc,=l 



1/2 



(23) 



where the number of replicas N^t for each value of the strong coupling is determined in 



rep 



the gaussian case by Eq. (14). In this case, we have taken as central value and uncertainty 
on a s those given in Eqs. (9 10) respectively. 



The results for the uncertainty obtained in this way are shown in Fig. [TTJ each normalized 
to the corresponding central cross section. They are compared to results obtained adding 
in quadrature the PDF uncertainties displayed in Fig. [7| and the a s uncertainties Eq. (20) 



displayed in Fig. 10, in turn obtained either with fixed PDFs, or by taking the PDF set that 
corresponds to each value of a s . Note that the range of a s variation for the MSTW08 curve 



with full correlation kept into account, given in Eq. (22), differs slightly from that, Eq. (10), 
used in all other cases. The effect of the correlation between a s and PDF uncertainties is 
indeed quite small: as one might expect, it is in fact smaller than the effect of the correlation 
between a s and PDF central values shown in Fig. 10 , and much smaller than the uncertainty 



on the uncertainty discussed in Sect. 4.2 
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6 The cross section: final results 



Our final results for the Higgs cross section are collected in Fig. [12] the cross sections are the 



same of Fig. [7j but now with the total PDF+a<, uncertainty. This is computed taking fully 
into account the correlation between a s , PDF central values and uncertainties for NNPDF 



and MSTW (as discussed in Sect. 5.3 and shown in Fig. 10). For CTEQ, which does not 



provide error sets for each value of a s , the PDF uncertainties (Fig. [7]) and a s uncertainties 



(Sect. 5.2 and Fig. 10) are added in quadrature; however, as seen in Sect. 5.3, this makes little 
difference. 

The main features of these final results are the following: 

• The total uncertainty on the cross section found by various groups are in reasonable 
agreement, especially if one recalls that they are typically affected by an uncertainty of 



order of about half their size, as seen in Sect. |4.2[ This follows from the fact that the 
total uncertainty is close to the sum in quadrature of the PDF and a s uncertainties, 
with the a s uncertainty essentially the same for all PDF sets, and PDF uncertainties 
on the input parton luminosity quite close to each other, 

The predictions obtained using the three given sets all agree within the total PDF+a s 
one-cr uncertainty. This is a consequence of the fact that the central values of the cross 
section computed using the same value of a s for all sets agree within PDF uncertainties, 
and that the spread of central values of a s used by the three groups Eq. ^ is essentially 
the same as the a s uncertainty under consideration. 



The results shown in Fig. 12 can be combined into a determination of the cross section 
and its combined PDF+a s uncertainty. We have seen that there is reasonably good agree- 
ment between the parton sets under investigation: apparent disagreement is only found if 
one compares results obtained with values of a s which differ more than the uncertainty on 
a s . However, in some cases (for example at the Tevatron for light Higgs) the agreement is 
marginal: the one-a uncertainty bands just about overlap. Ideally, this situation should be 
resolved by the PDF fitting groups by investigating the origin of the underlying imperfect 
agreement of parton luminosity. However, until this is done, a common determination of 
the cross section with a more conservative estimate of the uncertainty may be obtained by 
suitably inflating the PDF uncertainty. 

We have considered two different procedures which lead to such a common determination, 
based on the idea of combining a common a s uncertainty together with a PDF uncertainty 
suitably enlarged in order to keep into account the spread of PDF central values obtained 
using different PDF sets. 

Procedure A: This procedure is based on the observation that the change in PDFs and 
uncertainties when a s is varied by Aa s ~ 0.001 is small, as shown in Fig. [3j and 
implicitly demonstrated by the weak effect of the keeping it into account when evaluating 



a s and PDF uncertainties, Figs. 10 11 Therefore, we can obtain a prediction with a 
common value of a s for all groups by simply using the common intermediate value of 
a s = 0.119 in the computation of the hard cross section Eq. (|2|, and then using each 
group's PDFs and full PDF+a s uncertainties, despite the fact that strictly speaking they 
correspond to the slightly different values of a s listed in Eq. (J9J). Because all predictions 
are given at the same value of a s , their spread reflects differences in underlying PDFs. 
Hence, we can take as a conservative estimate of the one-cr total PDF+a s uncertainty 
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Figure 12: Cross sections for Higgs production from gluon-fusion the Tevatron (top), LHC 7 TeV 
(center) and LHC 14 TeV (bottom). All uncertainty bands are one-u combined PDF+a s uncertain- 



ties, as in Fig. [IT] (exact) for MSTW and NNPDF, and as in Fig. [10] for CTEQ, with the central 
value of a s of Eq. Q. The left column shows results normalized to the MSTW08 result, and the 
right column results normalized to each group's central result (relative uncertainty). 
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Figure 13: The CTEQ, MSTW and NNPDF curves obtained using the common value a s = 0.119 
in the cross section Eq. ([2]), but with the PDF sets and uncertainties corresponding to each group's 
value of a s Eq. All curves are normalized to the central MSTW curve obtained in this way. 
The common prediction can be taken as the envelope of these curves (Procedure A). The envelope 
curve shown (Procedure B) is instead the envelope of the CTEQ, MSTW and NNPDF predictions 
showed in Fig.[7j obtained including PDF uncertainties only, but with each group's value of a s used 
both in the PDFs and cross section, also normalized as the other curves. 



on this process the envelope of these predictions, i.e. the band between the highest 
and the lowest prediction. This procedure would be easiest to implement if all PDF 
groups were to provide PDFs with a common a s value and uncertainty. It is still viable 
provided the central a s values are not too different, and the a s uncertainty can be taken 
as the same or almost the same for all groups. 

Procedure B: This procedure is based on the observation that in fact the spread of central 
values of a s Eq. ^ is essentially the same as the width as the one-cx error band Eq. ( 10 ). 



Because the total uncertainty is well approximated by the combination of a s and PDF 
uncertainties, this suggests that we can simply substitute the a s uncertainty by the 
spread of values obtained with the three sets. Hence, a conservative estimate of the 
one-o" total PDF+a s is obtained by taking the envelope of the PDF-only uncertainty 
bands obtained using each of the three sets, each at its preferred value of a s Eq. ([9]), 
i.e. the envelope of the bands shown in Fig. [7} 



These two conservative estimates are shown in Fig. 13, where we display the uncertainty 
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bands obtained from the three MSTW, CTEQ and NNPDF sets whose envelope corresponds 
to the first method, as well as the envelope of the bands of Fig. [7j corresponding to the second 
method. The results turn out to be in near-perfect agreement, and we can take them as a 
conservative estimate of the PDF+a s uncertainty. Note that if any of the three sets were 
discarded, the prediction would change in a not insignificant way. 

Typical uncertainties are, for light Higgs, of order of 10% at the Tevatron and 5% at the 
LHC. Very large uncertainties are only found for very heavy Higgs at the Tevatron, which is 
sensitive to the poorly known large x gluon. As a central prediction one may take the midpoint 
between the upper and lower bands: in practice, this turns out to be extremely close to the 
MSTW08 prediction found adopting the previous method, i.e. using the MSTW08 PDFs but 
with a s = 0.119 in the matrix element. 

These results for the combined PDF+a s uncertainties should be relevant for Higgs searches 
both at the Tevatron and at the LHC. For instance, the latest combined Tevatron analysis on 



Higgs production via gluon-fusion 22 , which excludes a SM Higgs in the mass range 162-166 
GeV at 95% C.L., quotes a 11% systematic uncertainty from PDF uncertainties and higher 
order variations. It would be interesting to reassess the above exclusion limits if the combined 
PDF+a s uncertainties are estimated as discussed in this work. 
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7 Conclusions 



We have presented a systematic study of the impact of PDF and a s uncertainties in the total 
NLO cross section for the production of standard model Higgs in gluon-fusion. Whereas a 
full estimate of the uncertainty on this process would also require a discussion of other sources 
of uncertainty, such as electroweak corrections and the uncertainties related to higher order 
QCD corrections (NNLO, soft gluon resummation etc.) our investigation has focussed on 
PDF uncertainties, which are likely to be dominant for many or most LHC standard candles, 
and the a s uncertainty which is tangled with them. The process considered here is one for 
which these uncertainties are especially large, and thus it provides a useful test case. 
Our main findings can be summarized as follows: 

• Parton distributions are correlated to the value of a s in a way which is visible, but of 
moderate significance. In particular, if a s is varied within a reasonable range, not much 
larger than the current global uncertainty, uncertainties due to PDFs and the variation 
of a s can be considered to good approximation independent and the total uncertainty 
can be found adding them in quadrature. 

• The gluon luminosities determined from MSTW08, CTEQ6.6 and NNPDF1.2 agree to 
one a in the sense that their uncertainty bands always overlap (though just so for light 
Higgs at the Tevatron). As a consequence, the Higgs cross sections determined using 
these PDF sets agree provided only the same value of a s is used in the computation of 
the hard cross section. The spread of central values between sets is of the order of the 
PDF uncertainty. 

• The PDF uncertainties determined using these sets are in reasonable agreement and 
always differ by a factor less than two, while being affected by an uncertainty which is 
likely to be about half their size. The a s uncertainties are essentially independent of 
the PDF set, and provide more than half of the combined uncertainty. The combined 
uncertainties determined using the sets under investigation are in thus in good agreement 
with each other. 

• A conservative estimate of the total uncertainty can be obtained from the envelope of the 
PDF+a s uncertainties obtained from each set, all evaluated with a common central a s 
value. Equivalently, it can be obtained from the envelope of the PDF-only uncertainties 
of sets evaluated each at different value of a s , within a range of values which covers the 
accepted a s uncertainty. 

• A typical conservative PDF+a s uncertainty is, for light Higgs, of order 10% at the 
Tevatron and 5% at the LHC. This is at most a factor two larger than the PDF+a s 
uncertainty obtained using each individual parton set. Exclusion of any of the three 
sets considered here would lead to a total uncertainty which is rather closer to that of 
individual parton sets. 

Further improvements in accuracy could be obtained by accurate benchmarking and cross- 
checking of PDF determinations in order to isolate and understand the origin of existing 
disagreements. However, the overall agreement of existing sets appear to be satisfactory even 
for this worst-case scenario. 
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Note added: As this paper was being finalized, a study 23 of uncertainties on Higgs pro- 
duction has appeared. This paper presents detailed investigations of uncertainties which are 
not being discussed by us, specifically electroweak and scale (higher order QCD) uncertain- 
ties, while it addresses only marginally the issue on which we have concentrated, namely the 
interplay of a s and PDF uncertainties. 
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